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Increasing evidence suggests that single nucleotide polymorphisms (SNPs) associated 
with complex traits are more likely to be expression quantitative trait loci (eQTLs). 
Incorporating eQTL information hence has potential to increase power of genome-wide 
association studies (GWAS). In this paper, we propose using eQTL weights as prior 
information in SNP based association tests to improve test power while maintaining 
control of the family-wise error rate (FWER) or the false discovery rate (FDR). We apply 
the proposed methods to the analysis of a GWAS for childhood asthma consisting of 1296 
unrelated individuals with German ancestry. The results confirm that eQTLs are enriched 
for previously reported asthma SNPs. We also find that some SNPs are insignificant 
using procedures without eQTL weighting, but become significant using eQTL-weighted 
Bonferroni or Benjamini-Hochberg procedures, while controlling the same FWER or FDR 
level. Some of these SNPs have been reported by independent studies in recent literature. 
The results suggest that the eQTL-weighted procedures provide a promising approach for 
improving power of GWAS. We also report the results of our methods applied to the 
large-scale European GABRIEL consortium data. 

Keywords: asthma, family-wise error rate, false discovery rate, eQTL, genome-wide association study, weighted 
hypothesis test 



INTRODUCTION 

Asthma is a disorder characterized by inflamed mucosa of small 
airways of lung, causing wheezing and shortness of breath 
(Moffatt et al., 2010). Among the most common chronic dis- 
eases of childhood, asthma has been reported to affect more 
than 10% of children in many westernized societies (Cookson, 
2004). It is caused by a combination of genetic and environmental 
factors (Cookson, 2004; Moffatt et al., 2007), and several genome- 
wide association studies (GWAS) have been conducted to study 
the genetic basis underlying the complex disorder. More than 
50 single nucleotide polymorphisms (SNPs) have been reported 
to be associated with asthma, according to the GWAS catalog 
(www.genome.gov/gwastudies, accessed on January 15, 2013). 
Remarkably, the recent report (Moffatt et al., 2010) from the 
GABRIEL (A Multidisciplinary Study to Identify the Genetic and 
Environmental Causes of Asthma in the European Community) 
consortium identified several SNPs reaching genome-wide signif- 
icance through a large-scale meta-analysis. 

Prior biological information, often available in practice, has 
potential to increase power of GWAS. The common practice of 
GWAS, "agnostic" in some sense, assumes no prior informa- 
tion about any of the SNPs under investigation, meaning that all 
the SNPs have an equal likelihood of being causal. Some recent 



studies have taken advantage of information from linkage analysis 
(Roeder et al., 2006) and gene expression (Xiong et al., 2012) 
in genome-wide association scans. In genetic studies of etiol- 
ogy of asthma, it is of our particular interest to employ similar 
approaches and explore potentials of power gain in identifying 
asthma-associated SNPs by incorporating expression quantitative 
trait loci (eQTL) information. 

Catalogs of eQTLs in multiple tissues have been made publicly 
available, resulting from recent efforts of GWAS of gene expres- 
sions (Stranger et al., 2005, 2007, 2012; Dixon et al, 2007; Dimas 
et al., 2009; Yang et al., 2010). eQTLs provide insight into biol- 
ogy of transcription regulation. It has been shown that eQTLs 
are enriched for SNPs associated with complex diseases and traits 
using GWAS (Cookson et al., 2009; Nicolae et al., 2010). eQTL 
results can be used to provide functional interpretation for find- 
ings from GWAS (Moffatt et al., 2007; Heid et al., 2010; Hsu 
et al., 2010; Lango Allen et al, 2010; Speliotes et al., 2010; Chu 
et al., 2011; Wu et al., 2012) and prioritize genes in an associa- 
tion region for carrying out functional experiments using animal 
models (Teslovich et al., 2010). Focusing on eQTLs may also 
be useful to identify genetic pathways associated with the risk 
of complex diseases and traits, such as basal cell carcinoma in 
a skin cancer GWAS (Zhang et al., 2012) and type 2 diabetes 
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(Zhong et al, 2010). Other results show that many cis eQTLs are 
shared across tissues (Ding et al, 2010) and that a comprehensive 
eQTL catalog in one tissue might be used to increase the power of 
capturing relevant transcripts for other diseases (including those 
that are only weakly or incidentally expressed in tissues where 
eQTL information was collected). 

As single-SNP analysis still remains the most popular in 
GWAS, we focus on those methods designed for this type of 
analysis. Single-SNP analysis tests one SNP at a time for associ- 
ation by scanning across the whole genome, and hence involves a 
large number of hypotheses. To correct for multiple comparisons, 
statistical methods have been proposed and applied to control 
for the family-wise error rate (FWER) (Bonferroni, 1936; Holm, 
1979) or the false discovery rate (FDR) (Benjamini and Hochberg, 
1995; Storey and Tibshirani, 2003). Recent advances in statistical 
methodology make it possible to incorporate prior information 
through weighted hypothesis testing. In several of such methods 
(Genovese et al., 2006; Roeder et al., 2006, 2007), hypotheses are 
up-weighted or down-weighted based on prior likelihood of asso- 
ciation with the trait of interest. While keeping the FWER or FDR 
under control, the procedures can improve power with informa- 
tive weights and suffer small loss in power with uninformative 
weights (Genovese et al., 2006; Roeder and Wasserman, 2009). 
This feature is appealing as compared to prescreening SNPs based 
on prior information (e.g., to consider only eQTLs for associa- 
tion testing). In this paper, we propose to use eQTLs as prior 
information, and apply these weighted hypothesis testing meth- 
ods to reanalyze the MAGICS (Multicentre Asthma Genetics in 
Childhood Study) data of asthma GWAS (Moffatt et al., 2007) as 
well as the GABRIEL meta-study of asthma (Moffatt et al, 2010). 

RESULTS 

PUBLISHED ASTHMA ASSOCIATIONS ARE ENRICHED WITH eQTLs 

We extracted published asthma associations from the GWAS 
catalog maintained by the National Human Genome Research 
Institute. As of January 15, 2013, 52 distinct reference SNPs in or 
near more than 40 genes have been reported to be associated with 
asthma (Table Al). According to the eQTL database (described in 
Materials and Methods), 20 of these 52 SNPs (38.5%) are eQTLs. 
Using the proxy SNP search tool SNAP (Johnson et al., 2008), we 
then obtained an extended list of 506 SNPs that were either in the 
GWAS catalog or in strong linkage disequilibrium (LD) with the 
52 SNPs (r 2 > 0.8). We called all these 506 SNPs the extended set 
of asthma-associated SNPs. 

We calculated an eQTL enrichment p-value (Hosack et al., 
2003) using the MAGICS data. There are 300,821 SNPs that 
passed quality control in the MAGICS data. Among these SNPs, 
29 SNPs are in the GWAS catalog, and 64 SNPs are among 
the 506 extended asthma-associated SNPs defined previously. To 
account for the LD between SNPs in the calculation of enrich- 
ment p-value, we conducted LD pruning with the r 2 threshold 
of 0.8 on the 300,821 SNPs. This resulted in 251,826 SNPs and 
38 of them are extended asthma-associated SNPs according to 
the GWAS catalog. According to the eQTL database, 22,922 SNPs 
(9.1% of 251,826 SNPs) are eQTLs, and 13 asthma associated 
SNPs (34.2% of 38 SNPs) are eQTLs. The corresponding enrich- 
ment p-value is 6.78 x 10~ 5 , suggesting the asthma associations 



are enriched with eQTLs in the MAGICS data. Note that other 
analyses considered all the SNPs rather than the pruned set of 
SNPs. 

These results are in line with the previous findings (Nicolae 
et al., 2010), which studied the eQTLs in lymphoblastoid cell lines 
(LCL) from the HapMap samples and the GWAS catalog. Their 
results suggest that SNPs associated with complex traits are more 
likely to be eQTLs. 

WEIGHTS USING eQTL INFORMATION 

We calculated two kinds of weights using eQTL information for 
the MAGICS data. All the 300,821 SNPs passing quality control 
were considered. First, we defined a SNP as an eQTL SNP if 
it was labeled as an eQTL in the eQTL database (see details in 
Materials and Methods). There are 31,781 cis eQTL SNPs (10.6% 
of 300,821 SNPs) according to the definition, and for each of 
them, we retrieved an eQTL p-value p e QTL- Next, we considered 
two choices of weights, the general weight and the binary weight. 

The general weight is w g = ^j— log 10 p £ q TL for an eQTL SNP, 

and Wg = 1 otherwise. The binary weight takes only two possible 
values, Wf, = 3.70 for any eQTL SNP and wy = 0.68 otherwise. 
The two values of the binary weight were chosen to maximize 
the minimum power while keeping at least 10.6% (also the 
percentage of eQTL SNPs) of all the hypotheses with a power 
of 60%. The parameters for calculating the binary weight are 
6 = 0.106, a = 0.05, and f$ = 0.4 (see details in Materials and 
Methods). Last, both weights were normalized to have the mean 
equal to 1 which is necessary for the weighted hypothesis testing 
methods to maintain the correct FWER or FDR (Genovese et al., 
2006). After normalization, the general weight w g has a mean of 
2.44 and a median of 2.21 among eQTL SNPs, while the binary 
weight w b is 3.70 for all eQTL SNPs (Figure 1). 

WEIGHTED HYPOTHESIS TESTING 

We applied the weighted hypothesis testing methods (Genovese 
et al., 2006; Roeder and Wasserman, 2009) using the general 
weight w g and the binary weight Wb to the MAGICS data. For 
each of the 300,821 SNPs, we calculated the trait association p- 
value, p, from the single-SNP association test on the phenotypes 
of asthma status, as well as the weighted p-values Qg = p/w g and 
Q b = p/wy. Multiple testing adjustments were done for both the 
original p-values (p) and the weighted p-values (Q g and Q/,). 
Bonferroni (1936) and Holm's (1979) methods were considered 
to control for FWER, and Benjamini and Hochberg's (1995) 
method was used to control for FDR. 

We first ranked the SNPs using their p-values in the ascending 
order, and compared the ranks based on the weighted p-values 
with those based on the original p-values (Figure 2). Since only 
10.6% SNPs are eQTL SNPs according to the eQTL database, and 
hypotheses for eQTL SNPs are up-weighted, eQTLs generally have 
higher ranks after weighting, and non-eQTLs' ranks are lower but 
the magnitude of changes is small. This is true for both the general 
and binary weights. When restricting to the 29 asthma-associated 
SNPs reported in the GWAS catalog, we also observed similar 
behaviors, suggesting that weighting hypotheses may improve 
power using informative weights, and sacrifice a little power using 
uninformative weights (Roeder and Wasserman, 2009). 
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FIGURE 1 | Weights used in the MAGICS analysis. Each weight on the square root of — log 10 p e QTL where p e QjL is the eQTL p-value; (B) the 
corresponds to a SNP and a hypothesis. The weights have been normalized weights take only two possible values, which are decided using the method 
to have mean 1 and shown in the ascending order. (A) The weights are based described in Materials and Methods. 



We then looked at the Q-Q plots of the original and weighted 
p-values. For p-values greater than 0.0001, the Q-Q curves 
(Figure 3) are similar between the original and weighted p-values, 
regardless of the weights used. For those p-values less than 0.0001, 
some weighted p-values are smaller than original ones, and the 
difference is larger using the binary weight. We also observed that 
3 asthma-associated SNPs in the GWAS catalog are among the top 
SNPs with originalp-values less than 10~ 6 . The weighted p-values 
for all the 3 asthma-associated SNPs are smaller than original 
ones. 

Next, we applied the methods to control for the FWER. An 
effective ratio of 0.791 (Li et al., 2012) was used to calculate the 
effective number of SNPs (300,821 x 0.791). Controlling for an 
FWER level of 0.05, we obtained significant SNPs using both 
the original and weighted p-values. Both Bonferroni and Holm's 
methods gave the same results, and both weights (binary and gen- 
eral) also gave the same results (Tables A2, A3). The unweighted 
hypothesis testing claimed 6 SNPs to be significant, all on chro- 
mosome 17, including 2 asthma-associated SNPs (rs3894194 
with GSDMA, and rs7216389 with ORDML3) that have been 
reported previously (Moffatt et al., 2007, 2010). After applying 
the weighted hypothesis testing, we obtained 9 significant SNPs 
including all the 6 SNPs identified by the unweighted method, 
although the ranks are not exactly the same. The 3 SNPs addi- 
tionally identified by eQTL weighting were rs3902025, rs4795405, 
and rs2305480. The SNP rs2305480, a missense SNP in the gene 
GSDMB, was not reported in the previous GWAS study (Moffatt 
et al., 2007) but has been reported as an asthma-associated SNP 
in a later larger scale study by the GABRIEL consortium (includ- 
ing the MAGIC data, Moffatt et al, 2010) and was found to be 
strongly interacting with exposure to tobacco smoke in early life 
(Bouzigon et al, 2008). We found that rs2305480 is actually in 
LD (r 2 = 0.702, D' = 0.926) with rs7216389 that was identified 
by the unweighted methods, suggesting that rs2305480 may not 
represent a new association. Using a stringent r 2 threshold of 0.4, 
we found that the other two SNPs, rs3902025 and rs4795405, are 
also in LD with at least one SNP identified by the unweighted 



methods. So there is no new association identified by the weighted 
methods in this particular analysis. 

Besides controlling for FWER, we also used Benjamini and 
Hochberg's (BH) procedure (Benjamini and Hochberg, 1995) to 
control for a FDR level of 0.05 (Table A4). Based on the original 
p-values without weighting, the BH procedure gave 11 positive 
results (SNPs). The weighted BH procedures based on the general 
weight and the binary weight resulted in 7 and 8 additional posi- 
tive results (SNPs), respectively. Using a stringent r 2 threshold of 
0.4, we found that 5 SNPs (Table 1) are not in LD with any of the 
SNPs identified without weighting. Although none of the 5 SNPs 
are, or in LD with, any asthma-associated SNPs according to the 
GWAS catalog, there are some SNPs that seem interesting. Some 
of the SNPs are in or close to the genes PGAP3 and STARD3 on 
chromosome 17, and interestingly, rs2941504 has been reported 
in a recent independent study (Anantharaman et al., 2011) to 
be associated with asthma, although it does not meet the cri- 
teria for inclusion in the GWAS catalog. This suggests that the 
reanalysis using the eQTL weighting approaches is promising and 
potentially useful. 

REANALYSIS OF THE GABRIEL DATA 

As another application, we reanalyzed the GABRIEL data using 
the eQTL weighted approaches. Since only the p-values are nec- 
essary for the use of eQTL weighting, we took the p-values of the 
meta analysis of 37 studies that were calculated based on imputed 
data. In total, there are 2,473,850 SNPs and their p-values avail- 
able in the GABRIEL study, which include 267,350 out of 268,204 
eQTL SNPs in the eQTL database. The weights based on eQTL 
information were calculated in the similar way to the MAGICS 
data analysis. 

We applied Bonferroni and Holm's methods with an FWER 
level of 0.05, as well as the BH procedure with an FDR level of 
0.05. An effective ratio of 0.30 (Li et al., 2012) was used to cal- 
culate the effective number of SNPs (2,473,850 x 0.30). After 
obtaining the lists of significant SNPs using different methods, 
we report any SNPs identified by eQTL weighting that are not in 
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FIGURE 2 | Rankings of the SNPs based on original p-values and 
weighted p-values in the MAGICS analysis. (A) Original ranks of 
eQTLs compared to their new ranks based on the general weight; 
(B) original ranks of non-eQTLs compared to their new ranks based 
on the general weight; (C) original ranks of eQTLs compared to their 



new ranks based on the binary weight; (D) original ranks of 
non-eQTLs compared to their new ranks based on the binary weight. 
The black circles represent the reported asthma-associated SNPs in the 
GWAS catalog, and the gray circles represent the rest of the SNPs in 
the data. 



LD with any SNPs identified by unweighted methods using an r 2 
threshold of 0.4. Such SNPs may be informative and suggest new 
associations. Tables 2 and 3 show the SNPs that were identified 
based on the general weight and the binary weight, respectively. 

SIZE SIMULATIONS ON FWER 

We conducted simulations using 5000 permutations based on the 
MAGICS data, and calculated the percentage of having at least 
one false positive claimed by Bonferroni and Holm's methods 
(ot = 0.05, with an effective ratio of 0.791). In fact, any SNPs 
claimed significant using the two methods would be a false posi- 
tive. The calculated percentages (Table 4) provide estimates of the 
FWER. Bonferroni and Holm's methods give the same results. 
The results suggest that, under the null hypothesis, the FWER 



level is controlled for the methods based on both the original 
(unweighted) and the weighted p-values. The simulations con- 
firm the validity of the weighted hypothesis method (Genovese 
et al, 2006). 

DISCUSSION 

It is of substantial interest to enhance the power for identify- 
ing associations in the era of post-GWAS. Besides meta-analysis 
that has been proved successful in power gain (Moffatt et al., 
2010), incorporating prior information has also received increas- 
ing attention. Such information can be obtained from various 
sources and levels, such as linkage analysis (Roeder et al., 2006), 
gene expression (Yang et al., 2010), and annotation information 
of variants (Adzhubei et al., 2010), genes (Saccone et al., 2007), 
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FIGURE 3 | Q-Q plots of original p-values and weighted p-values in the MAGICS analysis. The weighted p-values are based on (A) the general weight, or 
(B) the binary weight. The reported asthma-associated SNPs in the GWAS catalog are shown in circles. 



Table 1 | Additional significant SNPs or positive results identified by eQTL weighting methods after accounting for linkage disequilibrium in 
the MAGICS analysis. 



Chr 


SNP 


Gene 


p- value 


Q g 


Qb 


Method 


17 


rs1 877031 


STARD3 


4.32 x 10" 6 


1.54 x 10" 6 


1.17 x 10" 6 


BH 


17 


rs931992 


TCAR STARD3 


5.61 x 10" 6 


1.96 x 10" 6 


1.52 x 10" 6 


BH 


17 


rs1 565922 


PGAP3 


7.28 x 10" 6 


2.55 x 10" 6 


1.97 x 10" 6 


BH 


17 


rs2941504 


PGAP3 


7.64 x 10" 6 


2.68 x 10" 6 


2.06 x 10" 6 


BH 


10 


rs11191325 


SUFU 


8.49 x 10" 6 


3.37 x 10" 6 


2.29 x 10" 6 


BH 



Q g and Qt, are the eQTL weighted p-values based on the general weight and the binary weight, respectively. Using both weights result in the same set of SNPs. 



Table 2 | Additional significant SNPs or positive results identified by eQTL weighting methods after accounting for linkage disequilibrium in 
the GABRIEL analysis. 



Chr 


SNP 


Gene 


p- value 


Q g 


Q b 


Method 


5 


rs244749 


MEF2C 


6.25 x 10" s 


2.51 x 10" 5 


1.54 x 10" s 


BH 


5 


rs10075941 




6.13 x 10" 5 


3.03 x 10" 5 


1.51 x 10" 5 


BH 


17 


rs7503195 




7.66 x 10" s 


3.07 x 10" s 


1.88 x 10" s 


BH 


17 


rs17637472 




6.13 x 10" 5 


3.56 x 10" s 


1.51 x 10" s 


BH 


9 


rs7047575 




6.87 x 10" s 


3.92 x 10" s 


1.69 x 10" 5 


BH 



The results are based on the general weight. 



and pathways (Wang et al., 2007). The so-called "agnostic" GWAS 
may benefit from incorporating useful prior information. In our 
study of asthma, gene expression information is of particular 
interest, as a recent study (Moffatt et al., 2007) identified several 
eQTLs associated with asthma. 

In the reanalysis of the MAGICS data (Moffatt et al, 2007), we 
applied recently developed statistical methods that can improve 
power by weighting hypothesis (Genovese et al., 2006; Roeder 
and Wasserman, 2009). Using eQTL information obtained from 
an independent dataset, we employed weighted procedures that 



up-weighted eQTL SNPs and down-weighted non-eQTL SNPs 
while controlling for the FWER or the FDR. It has been proved 
(Genovese et al., 2006) that any set of nonnegative weights 
can guarantee substantial power gain given informative weights 
and little power loss for uninformative weights. The property 
implies that the weighted procedures are robust to informa- 
tiveness of weights and to the uneven coverage of genes and 
expression targets on the genome. We took advantage of this 
robustness and applied the procedures to an asthma study. We 
found additional SNPs that were significantly associated with 
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Table 3 | Additional significant SNPs or positive results identified by eQTL weighting methods after accounting for linkage disequilibrium in 
the GABRIEL analysis. 



Chr 


SNP 


Gene 


p-value 


Qg 


Ob 




Method 


5 


rs736801 




2.15 x 10" 7 


8.81 x 10" 8 


5.29 x 


10- 8 


Bonferroni, Holm 


6 


rs2596450 


HCG26 


2.72 x 10" 7 


1.16 x 10" 7 


6.69 x 


10- 8 


Bonferroni, Holm 


5 


rs10075941 




6.13 x 10" 5 


3.03 x 10" 5 


1.51 x 


10" 5 


BH 


17 


rs17637472 




6.13 x 10" 5 


3.56 x 10" 5 


1.51 x 


10" 5 


BH 


5 


rs244749 


MEF2C 


6.25 x 10" 5 


2.51 x 10" 5 


1.54 x 


10" 5 


BH 


9 


rs7047575 




6.87 x 10" 5 


3.92 x 10" s 


1.69 x 


10" 5 


BH 


17 


rs7503195 




7.66 x 10" 5 


3.07 x 10" 5 


1.88 x 


10" 5 


BH 


6 


rs9273363 


HLA_DQB1 


8.38 x 10" 5 


4.41 x 10" 5 


2.06 x 


10" 5 


BH 


5 


rs4351182 




9.98 x 10" 5 


4.91 x 10" 5 


2.45 x 


10" 5 


BH 


2 


rs13391794 




1.01 x 10" 4 


5.62 x 10" 5 


2.49 x 


10" 5 


BH 


5 


rs10044342 


MEF2C 


1.10 x 10" 4 


4.59 x 10" 5 


2.71 x 


10" 5 


BH 


2 


rs6751196 




1.14 x 10" 4 


6.22 x 10" 5 


2.80 x 


10" 5 


BH 


6 


rs176095 


PBX2, GPSM3 


1.25 x 10" 4 


5.12 x 10" 5 


3.08 x 


10" 5 


BH 


2 


rs2675073 




1.34 x 10" 4 


5.40 x 10" s 


3.29 x 


10" 5 


BH 


z 


r s i y i ooz i 




1 v in-4 
1 .^D X IU 


7 QO \, m — 5 
/ .aZ. X I U 


3.60 x 


10" 5 


d n 


2 


rs10497621 


NUP35 


1.59 x 10" 4 


7.55 x 10" 5 


3.91 x 


10" 5 


BH 


5 


rs244750 


MEF2C 


1.70 x 10" 4 


7.08 x 10" 5 


4.17 x 


10" 5 


BH 


6 


rs9366689 


POM121L2 


1.72 x 10" 4 


8.73 x 10" 5 


4.23 x 


10" 5 


BH 


6 


rs7775759 




1.81 x 10" 4 


6.43 x 10" 5 


4.46 x 


10" s 


BH 


6 


rs7741091 




1.81 x 10" 4 


6.43 x 10" 5 


4.46 x 


10" s 


BH 


8 


rs6601649 




1.90 x 10" 4 


9.68 x 10" s 


4.68 x 


10" s 


BH 


6 


rs204994 




1.94 x 10" 4 


8.02 x 10" s 


4.76 x 10" 5 


BH 



The results are based on the binary weight. 



Table 4 | Family-wise error rate estimates in 5000 permutations. 



P-value 


FWER 


Original 


0.0454 


Weighted by general weight 


0.0458 


Weighted by binary weight 


0.0460 



Both Bonferroni and Holm's methods gave the same results in the same 
scenarios. 



asthma according to the weighting hypothesis methods. Some of 
them were interesting after we accounted for LD and compared 
them to literature. Our analysis was the first application of 
this approach to asthma GWAS studies, and the results suc- 
cessfully illustrated the use of eQTL weighting in the context 
of asthma studies. As another application, we also reanalyzed 
the GABRIEL meta-analysis p-values and reported corresponding 
results. 

It is noted that the weighted procedures can utilize eQTL 
information from a reference database. Multiple choices of eQTL 
databases have already been made available (e.g., Yang et al., 
2010; Liang et al., 2013), and future efforts may provide even 
better reference of eQTL information. For example, the eQTL 
information considered in our reanalysis was obtained through 
a single platform (Affymetrix HG-U133 Plus 2.0), and better 
coverage of gene expression profiling may be achieved through 
RNA-Seq technologies or by combining information from various 
platforms. 



Besides the weighted procedures, an alternative method of 
using eQTL information is to simply test association between 
eQTL SNPs and the trait of interest. Such a method is not rec- 
ommended in GWAS as it excludes non-eQTL completely and 
relies on the prior information too heavily. By contrast, weighted 
procedures make it possible to consider eQTLs and non-eQTLs 
simultaneously. More importantly, they can possibly increase 
power if the prior information is useful and are able to maintain 
the type I error under the null. 

Applying the weighted procedures in our reanalysis only 
requires p-value of eQTL SNPs. This flexibility means that such 
analyses can be applied to any existing GWAS data, even if they 
do not have accompanying gene expression data. Although gene 
expression may have tissue-specific patterns, a substantial frac- 
tion of eQTLs may be shared across tissues (Ding et al., 2010). 
Hence eQTLs developed from tissues that are not directly relevant 
to the outcome of interest, such as those from publicly available 
eQTL databases based on LCL, can be used to improve power on 
GWAS. It is possible that using eQTL information from relevant 
tissues may result in even more power gain, if such information is 
available. 

Besides the particular weighting hypothesis method (Genovese 
et al., 2006) we adopted, Bayesian methods are potentially alter- 
native strategies to incorporate eQTL information. The use of 
Bayes factors has been applied to genetic association studies 
(Wellcome Trust Case Control Consortium, 2007; Stephens and 
Balding, 2009). In single SNP analysis, a prior is assumed for 
each SNP effect [e.g., N(0, 0.2 2 ) under a model of association in 
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Wellcome Trust Case Control Consortium (2007)]. eQTL infor- 
mation can be naturally incorporated into the prior, although 
it may be challenging to choose a realistic yet tractable alter- 
native model and to assess error rates (Hoggart et al., 2008), 
especially with the eQTL weight. One possible choice is through 
modifying the variance of the prior, for example assuming a 
prior N(0.2 2 w), where w is a weight of eQTL signal. Another 
possible choice is to keep the variance the same and increase 
the probability of association for eQTLs a priori. It is of inter- 
est in future research to explore these possibilities and con- 
sider the extension of Bayesian methods to incorporate eQTL 
information. 

In our analysis, we took into account the LD between SNPs 
by considering the effective number of SNPs (Li et al., 2012). 
As an alternative, testing SNP sets for association has potential 
of improving power and reducing the correlation between tests. 
Since the focus of this paper is to demonstrate the use of eQTL 
information in association testing, we will consider the weighted 
correlated hypothesis in future research. 

Two choices of weights were applied in our analysis includ- 
ing a binary weight and a weight using strength of eQTLs, and 
the results using the two weights were similar in our analysis. 
Theoretical results exist (Roeder and Wasserman, 2009) for the 
optimal binary weight, which provide guidance in choosing the 
values of the weight. The weight taking advantage of the eQTL 
strength may possibly provide more useful information, and what 
is the best choice of weights is still under research. 

Through an application to an asthma GWAS, we demon- 
strated the usefulness of eQTL weights in GWAS. Although 
results may vary depending on the traits of interest and the 
underlying biological mechanism, the potentials of increas- 
ing power and little investment required for reanalysis make 
the eQTL-weighted procedures desirable for reanalysis of exist- 
ing GWAS data and useful for design and analysis of future 
studies. 

MATERIALS AND METHODS 

THE MAGICS ASTHMA GWAS SAMPLES AND DATA 

The MAGICS (Multicentre Asthma Genetics in Childhood Study) 
study data (Moffatt et al, 2007), part of the GABRIEL con- 
sortium, were reanalyzed by incorporating eQTL information. 
Quality control procedures were conducted similarly to a pub- 
lished protocol (Anderson et al., 2010). Individuals with missing 
phenotypes, elevated missing rates (>5%), or outlying heterozy- 
gosity rate were removed. Markers with an excessive missing rate 
(>5%), low MAF (<5%), or failing in the HWE test (p-value < 
10~ 5 ) were all excluded as well. The remaining dataset contains 
1296 individuals (647 affected and 649 unaffected) genotyped 
across 300,821 SNPs. 

To account for possible divergent ancestry and population 
stratification, principal component analysis (PCA) was con- 
ducted using EIGENSOFT 4.2 (Patterson et al, 2006; Price et al., 
2006). The genotype data were pruned for LD prior to the PCA. 
The PCA result (Figure Al) suggests that no obvious stratifica- 
tion exists, and the signal of the first principal component is 
very weak. In the subsequent analysis, we still included the first 
principal component as a covariate. 



LD pruning was considered only in the calculation of enrich- 
ment p-value. It was conducted using PLINK (vl.07, downloaded 
from http://pngu.mgh.harvard.edu/purcell/plink/) (Purcell et al., 
2007). A moving window with a width of 50 SNPs and a step 
size of 5 SNPs was considered, and pairwise LDs were calcu- 
lated and pruned if r 2 > 0.8 (corresponding PLINK arguments: 
"-indep-pairwise 50 5 0.8"). 

THE GABRIEL META-ANALYSIS p-VALUES 

Association testing results, including SNP ID and p-values, were 
obtain from a reanalysis of the GABRIEL consortium data using 
imputed SNPs (Bouzigon et al., personal communication). The 
meta-analysis considered imputation of SNP genotypes using the 
HapMap 2 reference data for 37 studies, and calculated a meta- 
analysis p-value for each SNP using available data. Imputed SNPs 
were kept for analysis if their imputation scores (Rsq) were >0.5 
and if their minor allele frequencies were > 1%. In total there were 
2,473,850 SNPs that passed the quality control. Only the SNP ID 
and the p-values of these SNPs were obtained and used for the 
reanalysis described in this paper. 

EXPRESSION QUANTITATIVE TRAIT LOCI DATA 

An eQTL database (http://www.hsph.harvard.edu/liming-liang/ 
software/eqtl/) resulting from an independent dataset was used 
as prior information to be incorporated in the GWAS. The sam- 
ple contains 405 siblings from a panel of families of British 
descent (MRC-A) (Dixon et al., 2007). Global gene expression 
in LCLs was measured using Affymetrix HG-U133 Plus 2.0 
chips. All siblings were genotyped using the Illumina Sentrix 
HumanHap300 BeadChip (ILMN300K) and/or the Illumina 
Sentrix Human-1 Genotyping BeadChip (ILMN100K). The SNP 
genotype data were further imputed using the MaCH program, 
and each SNP was tested for association with probes in the 
gene expression data. Restricting to cis eQTLs ( 1 Mb region) 
and controlling for the FDR of 1%, there are 515,947 tests 
with logarithm of odds (LOD) scores greater than 3.172, corre- 
sponding to 268,204 unique SNPs. In case a SNP has multiple 
p-values reported for associations with different probes, the min- 
imal p-value was used for that SNP. These 268,204 SNPs are 
considered as eQTLs, and the database contains information 
of their physical positions, LOD scores, p-values, and resid- 
ing or nearby genes. Details of the database are described by 
Liang etal. (2013). 

GENETIC ASSOCIATION ANALYSIS 

Genetic association analysis of the MAGICS data was conducted 
in PLINK. Logistic regression was used to test for disease-trait 
SNP association while adjusted for gender and the first princi- 
pal component. Meta-analysis on GABRIEL data was carried out 
by combining association results from 37 studies using a ran- 
dom effect model, and all computations were done using Stata 
software. 

p-VALUE WEIGHTING METHODS 

Consider m hypotheses H\, . . . , H m and their test statis- 
tic p-values, Pi,...,P m . Suppose there are weights 
W\ , . . . , W m available for the m tests, respectively, satisfying 
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Wi > 0 and ££i W, = m. Define Q, = P;/W; and let Q ( i) < 
• • • 5 Q(m) be the sorted values. Let P(i), ...,P( m ) an d 
W(!), . . . , W (m ) be the values in the corresponding order. Q, is 
sometimes referred to a "weighted p-value" (e.g., Roeder and 
Wasserman, 2009), although it is not ap-value. 

The weighted Bonferroni procedure is to reject any hypothe- 
sis Hj (1 < j < m) that satisfies Qj < ct/m, where a is the desired 
level of FWER. Genovese et al. (2006) showed that this procedure 
controls FWER at level no greater than a. 

Holm's weighted procedure (1979) is carried out as fol- 
lows: given the desired a level of FWER, if Q(i) > a/m, no 
hypothesis is rejected; otherwise, find the largest j that satis- 
fies Q(,-) < a/ J2k=i for all i < j, and reject the hypotheses 
corresponding to the j smallest Q/'s. Genovese et al. (2006) 
also prove that this procedure can work for a general setting 
of weights. 

We also consider Benjamini and Hochberg's procedure (1995) 
for controlling FDR. Given the desired level a, find the largest j 
such that Q(j) < a ■ j/m, and reject the hypotheses correspond- 
ing to the j smallest Q/'s. Genovese et al. (2006) prove that this 
procedure controls FDR at level a. 

eQTL INFORMATION AS WEIGHTS 

The eQTL p-values were used to construct weights for the SNPs 
in the asthma GWAS reanalysis. We considered two kinds of 
weights, the binary weight wy and the general weight w g . The 
binary weight takes only two possible values that are predefined, 
denoted by w e QTL and w non _ e QTL- For m hypotheses, a binary 
weight is defined as wi, = (wj, i, . . . , w\,, m ) where wj,j = w c qxl 
if the jth SNP is an eQTL SNP, and wj, ; = w non _ e QTL if it is 
not an eQTL SNR Given the values of a, p, and €, the opti- 
mal values of w e QTL and w non _ e QTL were chosen (Roeder and 
Wasserman, 2009) to maximize the minimum power among all 
the hypotheses while having at least a fraction e with high power 
1 — p. Here a is either the level of FWER or FDR. We also con- 
sidered a general weight, where the weight w g = (w gt i, . . . , Wg >m ) 
has w g j = ,J — log la p e QTL if the ;'th SNP is an eQTL SNP with the 
eQTL p-value£eQTL> and w g , } = 1 otherwise. The particular form 
was intuitively chosen prior to the reanalysis of the GWAS data 
in consideration of avoiding up-weighting top eQTL SNPs too 
much. Both Wb and w g were then normalized such that the means 
equal to 1, i.e., wb = 1 and wg = 1. 

REPORTED ASSOCIATIONS IN THE GWAS CATALOG 

Asthma-associated SNPs and genes reported in publications 
were retrieved from the online catalog of published GWAS 



on January 15, 2013. The catalog limits the associations to 
those with p-values less than 1.0 x 10~ 5 and records only one 
SNP with a gene or region of high LD unless there was evi- 
dence of independent association. The reported associations were 
compared against the findings in the asthma GWAS data we 
reanalyzed. 

LINKAGE DISEQUILIBRIUM INFORMATION 

To account for LD between SNPs, LD information based on 
HapMap 2 was obtained. The SNAP proxy search tool (http:// 
www.broadinstitute.org/mpg/snap/ldsearch.php) was used to 
obtain the information, based on the HapMap 2 (rel22) reference 
and a distance limit of 500kb. 

SIZE SIMULATION USING THE ASTHMA GWAS DATA 

Besides analyzing the MAGICS asthma GWAS data, we also 
conducted size simulations by permuting the disease status 
in the data. Logistic regression was considered where the 
dependent variable was the disease status (affected or unaf- 
fected) and the independent variables included a single SNP 
effect, gender, and the first principal component. The regres- 
sion was applied to all the ~300,000 SNPs across the whole 
genome. Five thousand permutations were done by permut- 
ing the disease status among all the individuals, and then 
the model was refitted for each SNP. In the end of simula- 
tions, 5000 permutation p-values were obtained for each of 
the ~300,000 SNPs. 
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RESOURCES 

R functions have been made available for users' convenience 
to compute the eQTL weighted p-values in order to conduct 
weighted procedures such as Bonferroni, Holm's and Benjamini- 
Hochberg procedures. The functions can be found at http://www. 
hsph.harvard.edu/liming-liang/eqtl-weighted-gwas. 
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APPENDIX 
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FIGURE A1 | Principal components analysis for the MAGICS study. The scatter plots of (A) PC1 vs. PC2, and (B) PC2 vs. PC3 are shown, with each 
circle corresponding to an individual in the asthma GWAS. 
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Table A1 | Published asthma associated SNPs in the GWAS catalog as of January 15, 2013. 



Region 


Chr 


SNP 


Context 


Gene 


eQTL 


1 q21.3 
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rs41z9zb/ 


Intron 


IL6R 
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1 


rsll 01999 


Intron 


DVLJ /AM 
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z 


rsyoU/yoy 


IntergGnic 


II 1QR1 II 1D1 
IL Ion I ,IL I n I 




zq i z. i 
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I SO //II OU 
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rso / / 1 1 bo 


Intron 


// 1 OR 1 

IL 1 on 1 


v 
Y 


4q31.21 


4 


rs /bobbbU 


Intergenic 


LUL/zyb/b 




4q31.21 


4 


rsooObzob 


Intron 


UAb i 




5q31.3 


5 


rsbob/y lo 


Intergenic 


NUrlr 1 




5q31.1 


5 


,-11 -7yi rr rro "7 

rs 1 1 /4bbo/ 


Intron 


C5orf56 


v 
Y 


bqzz. I 


o 


rs I oo/zbo 


Intergenic 


TCI D 

1 oLr 




5q31.1 


5 


rs I zybbob 


Intron 


II 1 o 

IL lo 




5q31.1 


5 


rszU/ob4o 


Intron 


CI /~"7 1 A £7 

oLLzzAb 


v 
Y 


bqo 1. 1 


o 


rszz44U Iz 


Intron 


nAUbU 


v 
Y 


5q1 2.1 


5 


rs I boozbb 


Intron 


rUL4U 




bpz I.OZ 


R 
D 


rsyzbob I b 


Intergenic 


RTMI O Ul A DRA 

d 1 NLz, rlLA-UnA 




6q27 


6 


rsb4bbU4z 


Intergenic 


T 
1 




6p21.32 


6 


rsybuuyz / 
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n LA-UUA 


Y 


6p21.32 


6 


rsyz /bbyy 


Intergenic 


Ul A Hi/~)A 1 

nLA-UUAz 
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6p21.32 


6 


rs/ //bzzo 
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ui a nno i 
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RTMI 1 
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D 
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opz l.oz 


R 
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Intron 
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v 
Y 
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6 
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Y 
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Q 
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CA A a r\o 
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rs I /bzb4/z 
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Intergenic 
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Y 


17q21.1 


17 


rs3894194 


Missense 


GSDMA 


Y 


17q12 


17 


rs2305480 


Missense 


GSDMB 


Y 


17q12 


17 


rs72 16389 


Intron 


0RMDL3 


Y 


19q13.42 


19 


rs1 6984547 


Intron 


ZNF665 




20p13 


20 


rs4815617 


nearGene-5 


KIAA1271 




22q12.3 


22 


rs2284033 


Intron 


IL2RB 





"Y" in the eQTL column means the corresponding SNP is an eQTL SNP according to the eQTL database described in Materials and Methods. "nearGene-5" is an 
NCBI dbSNP function code, meaning that SNP is S to and 2kb away from a gene. 
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Table A2 | The significant SNPs identified by the unweighted method and the weighted methods (using general weight or binary weight) in 
the MAGICS analysis. 



Rank 




Unweighted method 






General weight 






Binary weight 






Chr 


SNP 


p- value 


Chr 


SNP 


Q g 


Chr 


SNP 




1 

i 


1 7 


rcQQQ/1 1 OA * 

I Soc5y4 i y*+ 


/imp no 


1 7 


roQQQ/1 1 OA * 

r s ocsy4 i y^ 


i. do t-uy 


1 7 


roQQQ/1 1 OA * 

r sooy4 i y^ 


1 not na 
i .uo t-uy 


2 


17 


rs8079416* 


6.86E-09 


17 


rs8079416* 


2.70E-09 


17 


rs8079416* 


1.85E-09 


3 


17 


rs4795408* 


2.40E-08 


17 


rs4795408* 


9.66E-09 


17 


rs4795408* 


6.49E-09 


4 


17 


rs3859192 


3.58E-08 


17 


rs2290400* 


1.75E-08 


17 


rs2290400* 


2.00E-08 


5 


17 


rs2290400* 


7.39E-08 


17 


rs72 16389* 


1.76E-08 


17 


rs7216389* 


2.02E-08 


6 


17 


rs7216389* 


7.48E-08 


17 


rs3859192 


4.31 E-08 


17 


rs3859192 


5.25E-08 


7 








17 


rs3902025* 


1.05E-07 


17 


rs3902025* 


7.44E-08 


8 








17 


rs4795405* 


1.07E-07 


17 


rs4795405* 


9.83E-08 


9 








17 


rs2305480* 


1.40E-07 


17 


rs2305480* 


1.54E-07 


Bonferroni's correction is used with an FWER level of 0.05. "eQTL. 


Table A3 | The significant SNPs identified by the unweighted method and the weighted methods (using general weight or binary weight) in 


the MAGICS analysis. 
















Rank 




Unweighted method 






General weight 






Binary weight 






Chr 


SNP 


p- value 


Chr 


SNP 


Q g 


Chr 


SNP 




1 


17 


rs3894194* 


4.01 E-09 


17 


rs3894194* 


1.68E-09 


17 


rs3894194* 


1.08E-09 


2 


17 


rs8079416* 


6.86E-09 


17 


rs8079416* 


2.70E-09 


17 


rs8079416* 


1.85E-09 


3 


17 


rs4795408* 


2.40E-08 


17 


rs4795408* 


9.66E-09 


17 


rs4795408* 


6.49E-09 


4 


17 


rs3859192 


3.58E-08 


17 


rs2290400* 


1.75E-08 


17 


rs2290400* 


2.00E-08 


5 


17 


rs2290400* 


7.39E-08 


17 


rs72 16389* 


1.76E-08 


17 


rs7216389* 


2.02E-08 


6 


17 


rs7216389* 


7.48E-08 


17 


rs3859192 


4.31 E-08 


17 


rs3859192 


5.25E-08 


7 








17 


rs3902025* 


1.05E-07 


17 


rs3902025* 


7.44E-08 


8 








17 


rs4795405* 


1.07E-07 


17 


rs4795405* 


9.83E-08 


9 








17 


rs2305480* 


1.40E-07 


17 


rs2305480* 


1.54E-07 



Holm s method is used with an FWER level of 0.05. *eQTL. 
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Table A4 | The positive results (SNPs) given by the unweighted method and the weighted methods (using general weight or binary weight) in 
the MAGICS analysis. 



Rank 




Unweighted method 




General weight 






Binary weight 




Chr 


SNP 


p-value 


Chr 


SNP 


a g 


Chr 


SNP 




i 
i 


1 7 
I / 
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4.01 E-09 


1 7 


rc OQQ/l 1 Q/l * 

r sooy4 1 y4 
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i.uot-uy 


Z 
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l S4 / yo4uts 
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A 

4 


1 7 
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3.58E-08 


1 7 


r szzyu4uu 


1 7RP no 
I. / Dl-Uo 


1 7 
I / 


rcoocm/i nn* 
r szzyu4uu 


o nnp no 

Z.UUt-Uo 


O 


1 7 


rs2290400* 


7.39E-08 


1 7 


r s / z i Dooy 


1 7RP no 
I. /OlZ-Uo 


1 7 
I / 


r s / z i Dooy 


o nop no 

Z.UZ t-Uo 


P. 

o 


1 7 
I / 


rs72 16389* 


7.48 E-08 


1 7 


i sodoy i yz 


4. J 1 t-Uo 


1 7 
I / 


i s oooy i yz 


k of;p no 
o.zut-uo 


/ 


I / 


rs3902025* 


2.75E-07 


1 7 
I / 


rs4/yo4UD 


1 ri7P n7 
I.U/ t-U/ 


1 7 
I / 


rsjyuzUzo 


7 /i /i p no 
/.44t-Uo 


Q 
O 


I / 


rs4795405* 


3.64E-07 


1 7 
I / 


rsoyuzuzo 


1 ncp n7 
I.Uot-U/ 


1 7 
I / 


rs4/yo4Uo 


Q QQP no 

y .oot-Uo 


Q 

y 


I / 


rs2305480* 


5.68E-07 


1 7 
I / 


roOQnfi/i on* 
rSzJUo4oU 


1 /inp n7 
l.4Ut-U/ 


1 7 
I / 


rrOQnc/i on* 
rSzoUo4oU 


i e.a p n7 


IU 


1 / 


rs1 1557467 


2.19E-06 


1 7 
I / 


rsoUo/ J /o 


/I CQC H7 

4. Do t-U/ 


1 7 
I / 


roQnC7Q7Q* 

rsoUo/ J /o 


c /i CP n7 
o.4ot-U/ 


1 1 


I / 


rs8067378* 


2.02E-06 


1 7 
I / 
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rsy JUoz / / 


7Q/I P ri7 
/.y4t-U/ 


1 7 
I / 


rrQQnQT77* 


o Qnp n7 

y.oUt-u/ 


12 








17 


rs1877031* 


1.54E-06 


17 


rs1 877031* 


1.17E-06 


13 








17 


rs907092* 


1.56E-06 


17 


rs931992* 


1.52E-06 


14 








17 


rs931992* 


1.96E-06 


17 


rs907092 


1.68E-06 


15 








17 


rs2941504* 


2.68E-06 


17 


rs2941504* 


2.06E-06 


16 








17 


rs1 565922* 


2.55E-06 


17 


rs1 565922* 


1.97E-06 


17 








17 


rs 11 557467 


2.64E-06 


10 


rs11191325* 


2.29E-06 


18 








10 


rs11191325* 


3.37E-06 


17 


rs1 1557467 


3.22E-06 


19 














17 


rs1007654* 


3.11 E-06 



The BH procedure is used with an FDR level of 0. 05. "eQTL. 
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